Maple 2024 Questions and Posts

These are Posts and Questions associated with the product, Maple 2024

Here is a minor nit regarding combining plots in Maple.

restart;

kernelopts(version);

`Maple 2024.2, X86 64 LINUX, Oct 29 2024, Build ID 1872373`

(1)

 

Plot the cardioids  1+`&+-`(cos*t) in polar coordinates:

pic1 := plot(1-cos(t), t=0..2*Pi, coords=polar);
pic2 := plot(1+cos(t), t=0..2*Pi, coords=polar);

 

 

 

When combining the plots, display() inherits the plotting range from the previous

plots, incorrecting taking it to be 0 .. 2*Pi
 

plots:-display(pic1, pic2, scaling=constrained);

 

To get the correct plot, we need to override the range determined by display():

plots:-display(pics, scaling=constrained, view=[-2..2,-2..2]);

 
 

 

Download mw.mw

Hi,

For a pedagogical purpose, I am trying to build an animation that illustrates the connection between:

  • the area under the curve of a probability density function, and

  • the cumulative distribution function of a random variable.

The idea is to let a bound x vary and show simultaneously:

  • the accumulated area under the density curve up to x;

  • the corresponding value of the CDF .

My goal is to obtain a visual representation that helps students clearly understand that the cumulative distribution function represents the accumulated area under the density function.

Do you have suggestions on how to organize this cleanly in Maple, possibly with synchronized animation or an interactive visualization?

Thank you very much for your help.

Q_CDF_PDF.mw

I am studying a nonlinear wave equation and trying to reproduce the energy balance method shown in a research paper. First, the original partial differential equation is reduced to an ordinary differential equation using a traveling wave transformation. After obtaining the reduced equation, the paper rewrites it in a form suitable for the energy balance method and derives the corresponding variational principle and Hamiltonian invariant. Then a trial periodic solution in cosine form is assumed. Using the Hamiltonian invariant and some initial conditions, the parameters of the trial function are determined and a periodic solution is obtained.

I would like to know how to implement this procedure in Maple. Specifically,  compute the Hamiltonian invariant from the equation, substitute the cosine trial function, and determine the unknown parameter in the trial solution using the energy balance method. I will attach images from the paper that show the derivation steps I am trying to reproduce. Any guidance on how to perform these symbolic steps in Maple would be very helpful.

f-s.mw

Hi,

For a pedagogical purpose, I am trying to illustrate the orthogonal projection H of a point A on to a plane P1

I constructed the line l1​ passing through A and perpendicular to the plane P1
However, in the graphical visualization, the line does not appear to be perpendicular to the plane. Visually, it gives the impression that the line is not orthogonal to the plane.

Do you have any idea what might cause this effect?

Thank you for your help.

Q_Espace.mw

I have two equations developed in different ways and I want to compare them.
The first is a development done using maple, while the second is an algebraic analysis.

I need to find out if both are identical and, therefore, if the subtraction would be equal to zero.

How should I proceed?

I had thought about simplifying the first and comparing it with the second; is that a good approach?

 

restart

with(plots)

H0 := -S1^2*eta1-S2^2*eta2-S1*gamma1-S2*gamma2

-S1^2*eta1-S2^2*eta2-S1*gamma1-S2*gamma2

(1)

NULL

Z0 := exp(-beta*H0)

exp(-beta*(-S1^2*eta1-S2^2*eta2-S1*gamma1-S2*gamma2))

(2)

Z0 := add(Z0, S1 = [-2, -1, 0, 1, 2])

Z0 := add(Z0, S2 = [-2, -1, 0, 1, 2])

G0 := -ln(Z0^((1/2)*N))/beta

G01 := -ln((2*exp(4*beta*eta1)*cosh(2*beta*gamma1)+2*exp(beta*eta1)*cosh(beta*gamma1)+1)^((1/2)*N)*(2*exp(4*beta*eta2)*cosh(2*beta*gamma2)+2*exp(beta*eta2)*cosh(beta*gamma2)+1)^((1/2)*N))/beta

-ln((2*exp(4*beta*eta1)*cosh(2*beta*gamma1)+2*exp(beta*eta1)*cosh(beta*gamma1)+1)^((1/2)*N)*(2*exp(4*beta*eta2)*cosh(2*beta*gamma2)+2*exp(beta*eta2)*cosh(beta*gamma2)+1)^((1/2)*N))/beta

(3)

Teste := G0-G01

NULL

Download TESTE_MAPLE.mw

my question is a little bit long but is not complicated, i want find thus  unkown but  realy i am don't know how apply on it by maple, i have two best paper which explain very well i just want to find thus dimensional Lie algebra which is be invariant or not satisfy condition or not which i have to used or which i have not to use it also the importan part how find them in paper 1 first , for equation fisher 
How find eq(29) which is i think is two -dimensional Lie algebra of equation, also the best part is reduction which by apply this we can change PDE to ode but i don't know how apply eq(31) or even find it yet  is related to eq(27-28)  and by replacing equation eq(34) we can get our ode i am just loking for the ode, for the eq(76) and eq(85) have same procedures,  i will mention the paper link too 

Lie.mw

paper-1

paper-2

Hi,

I am trying to reproduce a polar curve, as shown in the attached picture.
I have tested several polar expressions and different plotting options in Maple, but I can’t manage to obtain an exact match

If you have any ideas about a possible polar equation, a combination of functions, or plotting parameters (sampling, domain, polarplot options, etc.), I would be very grateful.

Thanks in advance for your help!

Q_Polar.mw

I am trying to calculate the stability of my solution using this integral test:

U = 1/2 * ∫ u(x,t)^2 dx

Then I want to compute the derivative with respect to epsilon and check if the result is greater than zero or less than zero.

I tried to do it in Maple but it is not working correctly.

Also, I want to substitute random numbers for the parameters to test stability. I prefer to use simple integer values like 1, 2, 3, etc., not decimal numbers. When I try random values, sometimes the result does not evaluate or gives complicated output.

S-test.mw

In the attached file, I want to calculate the integral Q1. Numerically, this is easy to do in Maple. For theoretical reasons, the exact result pi/e is known. However, a contradiction arises between command lines (4) and (5). Command (6) is also unsuccessful, as its exact result is unknown. What am I doing wrong?

restart

Q1 := int(sin(Pi*x)/(x^x*(1-x)^(1-x)), x = 0 .. 1)

int(sin(Pi*x)/(x^x*(1-x)^(1-x)), x = 0 .. 1)

(1)

evalf[100](Q1)

1.155727349790921717910093183312696299120851023164415820499706535327288631840916939440188434235673559

(2)

evalf[100](Pi/exp(1))

1.155727349790921717910093183312696299120851023164415820499706535327288631840916939440188434235673559

(3)

is(Q1 = Pi/exp(1))

false

(4)

identify(evalf[100](Q1))

Pi*exp(-1)

(5)

identify(.2340257795502385151002175791580229871350403567739388325733228478980460706709848394726222465477567339)

.2340257795502385151002175791580229871350403567739388325733228478980460706709848394726222465477567339

(6)

NULL

Download test.mw

is been a while i work on it but i can't figure out where is problem and even my solution is so far from it when numerically i tested , i got same outcome as the paper did maybe is  long but it is same and maybe have some typical different but they are same , the problem in here is that which when i substitue is not my answer i don't know where is my mistake ?

pde-te.mw

f18.mw

f19.mw

The plot in the attached file only works if the complete function expression is entered. If only the function name is entered, no plot appears. What am I doing wrong?test.mw

restart;

c(x)dsum(sin(10^k*x)/10^k, k = 1 .. 10^2)

plot(c(x), x = -(1/4)*Pi .. (1/4)*Pi)

 

``

c(x)dsum(sin(10^k*x)/10^k, k = 1 .. 10^2)

plot(sum(sin(10^k*x)/10^k, k = 1 .. 10^2), x = -(1/4)*Pi .. (1/4)*Pi)

 

 

NULL

Download test.mw

Hi,

I’m trying to transpose an existing animation that connects the unit circle to the graphs of cos⁡(θ)\ and sin⁡(θ)\ into a complex-numbers visualization, so that students can clearly see the link between

z=eiθ=cos⁡θ+isin⁡θ,arg⁡(z)=θ,∣z∣=1

and the corresponding real/imaginary components.

Goal: a dynamic view where a point z(θ) moves on the unit circle in the complex plane while (simultaneously)

  • the projections show ℜ(z)=cos⁡θ\  and ℑ(z)=sin⁡θ,

  •  the graphs of cos⁡θ and sin⁡θ are traced against θ\,

  • and/or the angle θ\ and argument are displayed in a clean, didactic way.

To better illustrate my objective, here is the link to the target animation I would like to transpose: 

Illustration

Thank you in advance for your insights and feedback.

Animation_Question.mw

Hello everyone
Dear experienced and expert friends
As a beginner, I would like to ask if any of my friends can guide me.
The following commands are related to Mathematica:

plots = Table[n = sValues[[i]];
   ParametricPlot[{1 - 2/n - 1.5/n^2 + (1.33 - 2/n) \[Gamma] - 
      0.0740741 (15 + 4*n) \[Gamma]^2, 
     12/n^2 + (16 \[Gamma])/n + (80 \[Gamma]^2)/9}, {\[Gamma], 0, 
     0.06}, PlotStyle -> colors[[i]], 
    PlotRange -> {{-10, 10}, {-10, 10}}], {i, Length[sValues]}];

Show[plots, Frame -> True, FrameLabel -> {"\!\(\*
StyleBox[SubscriptBox[\"n\", \"s\"],\nFontSize->16,\n\
FontColor->GrayLevel[0]]\)", "\!\(\*
StyleBox[\"r\",\nFontSize->16,\nFontColor->GrayLevel[0]]\)"}, 
 GridLinesStyle -> Black, PlotRange -> {{0.94, 1}, {0, 0.06}}, 
 PlotLegends -> 
  Placed[LineLegend[sValues, LegendLabel -> "s,w"], {0.5, 0.5}], 
 ImageSize -> 400]

I want to rewrite this process in Maple for my own functions.
I would be grateful if it is possible or if these commands are rewritten in a complete and executable form in Maple for me so that I can understand the working pattern. Or at least an equivalent command that can do this in Maple is introduced
Thank you all

Hi,
I’m taking the liberty of getting back to you once again regarding the animation of a procedure (here, for fractals). I’m able to use the Explore command correctly, but I’m having difficulty with animate.
Thank you very much for your help — it is truly invaluable.

Q_Fractales.mw

Hi,
I’d appreciate your insights on how to animate this sequence (family) of functions in order to obtain a rendering like the one shown below.


Any suggestions or best practices would be very welcome. Thank you!

AnimQ.mw

1 2 3 4 5 6 7 Last Page 1 of 47